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The Heavy Fermion state in UPd2Al3 may be approxi- 
mately described by a dual model where two of the three 
U- 5f electrons are in a localized state split by the crystalline 
electric field into two low lying singlets with a splitting energy 
A ~ 6 meV. The third 5f electron has itinerant character and 
forms the Heavy Electron bands. Inelastic neutron scattering 
and tunneling experiments suggest that magnetic excitons, 
the collective propagating crystal field excitations of the lo- 
calized 5f electrons, mediate superconducting (sc) pairing in 
UPd2Al3 . A theory for this novel mechanism is developed 
within a nonretarded approach. A model for the magnetic ex- 
citon bands is analyzed and compared with experiment. The 
sc pair potential which they mediate is derived and the gap 
equations are solved. It is shown that this mechanism favors 
an odd parity state which is nondegenerate due to the com- 
bined symmetry breaking by the crystalline electric field and 
the AF order parameter. A hybrid model including the spin 
fluctuation contribution to the pairing is also discussed. 
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I. INTRODUCTION 

Among the uranium- based Heavy Fermion (HF) 
superconductorsti (sc) unconventional behaviour is the 
rule rather than the exception. In the canonical inter- 
metallic U- compounds one observes orbitallv- degener- 
ate sc states with associated double transitiona aad two 
upper critical field curves as in the case of UPtstI, fur- 
thermore coexistence and competition of different un- 
conventinal sc states with a. SDW- like state in the 
Ui_a;,Tha;Bei3 was reportedo. In addition URu2Si2 
shows superconductivity embedded in a xegion with a 
yet undetermined hidden order parameterO. In all these 
cases, when antiferromagnetism (AF) is observed it has 
no well defined long range order and is characterized by a 
very small moment size of order 10~^/is or less. The s: 
nation is distinctly different in the compound UPd2Al3i 
Its normal state is characterized by a specific heat 7- 
value of 120 mJ/molc much smaller than observed in 
the previous compounds, it is therefore of only moder- 
ate HF character and is superconducting below Tc=1.8 
K. On the other hand it has long range AF order below 
TAr=14.3 K with much larger moments of almost atomic 
like size with /i — 0.85fj,B- This indicates that in addition 
to the itinerant electrons there must be nearly localized 
5f- electrons present. They result from the 5f^ configura- 



tion of the U*^ ionic species which is the dominating one 
in this compoundQ. This is supported by the tempera- 
ture dependent susceptibility which shows an extremely 
strong hexagonal ac- anisotropy with an almost constant 
Xc and a much larger XaiT) which is strongly temper- 
ature dependent with a typical behaviour known from 
crystalline electric field (CEF)- split localized 5f states 
with a singlet- (nonmagnetic) ground state and a first ex- 
cited state at an energy estimated as A ~6 meV. Knight 
shift analysis in the normal state obtained from /isr- ex- 
periments confirm the presence of localized 5f momentso. 
Since the CEF ground state is a singlet the AF order 
must be of the induced moment type, i.e. due to mixing 
with the first excited CEF state caused by the inter- site 
exchange. 

The most direct confirmation for this dual nature of 5f- 
electrons in UPd2Al3 is obtained from inelastic neutron 
scattering experiments which have been able to identify 
the collective propagating excitations that originate in 
CEF excitations of energy A and are broadened into a 
dispersive band usually termed 'magnetic excitons' due 
to the action of the intersite exchange. These modes 
extend up to 10 meV and for a wave vector along the 
hexagonal axis are quitpj^arp and well definedEI. More 
recently it was foundE2l"t3 that at the AF wave vector 
Q — (0, 0, tt) in reduced units (App. A) a strong in- 
teraction of these collective modes of localized moments 
with the heavy conduction electrons exists which leads 
to a resonance- like structure in the dynamical structure 
function due to a near degeneracy of the exciton mode en- 
ergy at Q and the sunerconducting gap. Before that tun- 
neling measurementsEj have already shown that strong 
coupling anomalies in the tunneling current exist in the 
same energy range where the magnetic excitons at Q were 
found. This has lead to the conclusion that they are in- 
deed the bosonic 'glue' which is responsible for the forma- 
tion of Cooper pairs in this compoundO. The argument 
is quite analogous to the strong coupling conventional 
electron-phonon superconductors where the phonon spec- 
trum known from inelastic neutron scattering leaves its 
imprint on the sc tunneling spectrum. UPd2Al3 is there- 
fore the first case of an unconventional HF superconduc- 
tor where a similar comparison has lead to the identifi- 
cation of the pairing mechanism, in this case the pairing 
is mediated by the exchange of magnetic excitons. 

In the present work this novel type of pairing mecha- 
nism is investigated in detail. It is different from both the 
electron- phonon mechanism in conventional supercon- 
ductors and from the spin fluctuation mechanism which 
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is commonly assumed to be present in the unconventional 
superconductors with strongly correlated electrons. The 
former is mediated by slightly damped real frequency 
bosons (phonons) which do not couple to the spin de- 
grees of freedom, thus allowing only for spin singlet pairs. 
The latter is mediated by strongly overdamped spin- fluc- 
tuations, i.e. bosons with purely imaginary frequency 
which couple to conduction electron spins in a rotation- 
ally invariant manner and in principle allow for both spin 
singlet and triplet pairing. The new pairing mechanism 
discussed here has distinctly different features: It is me- 
diated by magnetic excitons which are real frequency, 
propagating bosonic modes that couple to the conduction 
electron spin, however in a way which explicitly breaks 
spin rotational invariance due to the presence of the CEF 
splitting. Together with the effect of the AF order pa- 
rameter this will lead to a complete splitting of triplet 
states resulting in nondegenerate odd parity pair states. 

It is the aim of this work to study within a generic 
dual model of 5f electrons (Sect. II) the magnetic exciton 
bands (Sect. Ill) and the associated pair potential (Sect 
.IV) of this mechanism, incorporating both the localized 
CEF states which originate from the 5f^ configuration of 
the U''+ state of UPd2Al3 and the conduction electrons. 
Furthermore the gap equations and their explicit solu- 
tions together with their node line structure for a quasi 
one dimensional model are discussed explicitly; it is also 
shown how the presence of the AF background influences 
the gap structure and lifts the degeneracy (Sect .V). Fi- 
nally we discuss a hybrid model including also the usual 
spin fluctuation contribution (Sect. VI) and a summary 
and outlook is given in Section VII. 

II. BASIC HAMILTONIAN FOR THE DUAL 5F- 
MODEL 

The pronounced susceptibility anisotropy of UPd2Al3 
with an almost constant Xc and a strongly temperature 
dependent Xa(T) with a maximum at T=50 K is per- 
haps the most diregt evidence for the presence of lo- 
calised CEF- statesQ. In fact the Xa,c(T)- dependence 
was used to extract the hexagonal CEF- scheme in Ref.( 
1^ where it was conluded that U^+(5f^) ground state and 
first excited states are both singlets. However later it 
was founda that this cannot explain the temperature de- 
pendence of the staggered magnetisation and a singlet- 
doublet system was proposed instead. The magnetic ex- 
citations are very similar for both CEF systems except 
for a small difference in their temperature dependence. 
For simplicity we therefore use the singlet-singlet sys- 
tem with the ground state \g) (0 meV) and excited state 
|e) (A= 6 meV). The value of A is obtained from a fit 
to the experimental excitations (Sect. IV) and is close to 
the one obtained from the temperatuije variation of sus- 
ceptibility and AF order parameteiO. All higher CEF 
levels starting at ~ 10 meV will be neglected. The total 
Hamiltonian, including the localized part IIc_bf + ^//j 



the conduction band states (He) of the remaining itiner- 
ant 5f- electron and its interaction with localized states 
(He/) was introduced in Ref.( ^ as 

H = Hc + HcEF + Hff + Hcf 

- E ^//(uVl/, -2/o(ff-l)^5,/i 

Here e^:^ is the dispersion of itinerant electrons de- 
scribed by (Cjr^,Si) creation and spin operators. The 
most important Fermi surface (FS) sheet has the shape 
of a, sfehtly corrugated cylinder along the hexagonal c- 
axisllSo. It can be modeled by the expression 

Cfe^ = e_L (fcj_cr) - 2i II cos fc^ (2) 

Where eji_{k^a) is the dispersion perpendicular to the 
c-axis whose precise form is unimportant in the following, 
the much smaller dispersion || to c which is responsible 
for the corrugation of the FS cylinder is determined by 
an effective hopping energy ty. The localised 5f- elec- 
trons have a total angular momentum Ji (J=4). Within 
the singlet-singlet subspace they can be represented by a 
pseudospin Si (S=|) with the correspondence 

Jx = aSx, Jy = aSy, Jz = e(- - Sz) (3) 

This means that {e\Jx\g)— -i{e\Jy\g)— \a and {e\Jz\g) 
=0, i.e. inelastic singlet-singlet transitions can only be 
excited by the transverse operators ix.y The third term 
describes a possible superexchange J//= 4(g-l)^t^y/U 
between the localised 5f- electrons due to their remaining 
finite inter- site hopping t//. The last term finally rep- 
resents the exchange interaction which couples the itin- 
erant and localised 5f- electrons of the dual model. The 
total effective inter- site exchange between localized 5f- 
moments is then given by its Fourier transform as 

^J{q)=Tff{q)+ll{g-l)\'M 

x'o{q) = Xo{q)-j^J2^o{Q) (4) 
q 

Where the second (RKKY) part is in principle de- 
termined by the static conduction electron susceptibility 

Xo{q)- However in the following J{q) will be treated as 
an empirical quantity to be parametrized by fitting the 
actually observed magnetic excitations to the theoretical 
predictions based on the dual Hamiltonian of Eq.(|^). 

III. THE ANTIFERROMAGNETIC 
SINGLET-SINGLET SYSTEM 

The 'induced' magnetism of singlet-singlet CEF- sys- 
tems is well studied for Pr- compounds (for a review, see 
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Ref.( ^0|)). Since U^"'" has also a P electron configura- 
tion similar to Pr"^^ singlet-singlet systems may also exist 
in U-compounds with rather localised 5f- states such as 
UPd2Al3 . In a singlet-singlet system the ground state 
is nonmagnetic, i.e. {g\J\g) = 0. Nevertheless magnetic 
moments may appear spontaneously at T^r if the effec- 
tive exchange is strong enough to mix the excited state 
|e) into the ground state \g). The mixed state \g) will 
then have nonzero moment due to the nondiagonal ma- 
trix element a=2{e\Jx\g) ■ In this case the AF transition 
is preceded by the softening of a 'magnetic exciton' mode 
at the AF wave vector [Q = (0, 0, tt) in UPd2Al3 ). 



A. Exchange model and origin of induced AF order 

The magnetic exciton mode can be interpreted as orig- 
inating from local CEF- excitations \g) <-> |e) with energy 
A which propagate from site to site due to the effective 
exchange and this process leads to the magnetic exci- 
ton dispersion. It is therefore a collective excitation of 
localized 5f- CEF states and may be obtained from the 
second and third term in Eq.(|]) but taking the total effec- 
tive exchange of Eq.(^) instead. Within the two singlet 
subspace this leads to a Hamiltonian 




H, = {^/e)Y,Jl-\Y,J^^J^JJ, 



(5) 



where J is the uniaxial exchange tensor of Eq.(P) with 



ixx _ 



Its Fourier transform for the 



two AF sublattices A,B (=A, r) is given by (ie A) 

Jxriq) = ^Jtj exp{-iq{Ri - Rj)) = Jrxiq)* (6) 

As mentioned in Sect. (II) the exchange functions will 
be empirically parametrized to fit the magnetic exciton 
dispersion. The exchange couplings included are illus- 
trated in Fig. |l|. This leads to 



Joiq) = J AA{q) = JBsiq) = 

JN{q) = J AB[q) = JBAiq)* = 

[Jo + J2'-fi{q)]io{q) 

The structure functions ji^q) are given by 
70 (g) = 2cosg2 

1 %/3 
7i(g) = 2[cosqj; + cos(-g^ + ~^y^ + 



(7) 



,1 
cos(-fe 

72(17) = 2cos2g^ 



%/3 



(8) 



V3, 



73(17) = 2 [cos Viqy + cos— {qy + V^Qx) + 
cos ^{qy - Viqx)] 

where 70 (gz±(5^) = -70(92) andj2iqz±Qz) = l2{qz)- 
This leads to the important symmetries J d{Qz ^ Qz) = 
Joiqz) and JN{qz ± Qz) = -Jniqz) where q and 
(5=(0,0,7r) are given in units of aT^ and c^^ for x,y and 
z components, respectively. Each of the exchange pa- 
rameters J u (i^=0-4) in Eq.(^ corresponds to a given in- 
or out of plane neighbor shell in Fig. ^ is in general a 
uniaxial tensor. The tensor notation will now be sup- 
pressed for simplicity. In mean field (mf) approximation 
the Hamiltonian in Eq.(B[) reads: 



FIG. 1. Schematic view of conventional unit cell of 
UPd2Al3 with only U- atoms shown (a=5.35 A, c=4.185 A). 
Arrows indicate magnetic moments with FM- in plane order 
and AF stacking along c corresponding to Q=(0,0,^). The 
exchange pairs used in the model of Sect. 3 are also indicated. 
J3 and J4 describe the coupling to n.n.n. pairs along c and in 
the ab- plane respectively. The values of Ji(i=0-4) are given 
in the caption of Fig. 4. 



[{A/e)JS)-KJx{i)] 



he = [Jaa{0) - Jab (0)] (J) A = J{Q){J)x (9) 
Je = J{Q) = 2[3Ji + 3J4 + J3 - Jo - 6J2] 

Here h^— he{A)= -hf,{B) =heX is the staggered molec- 
ular field and (J) = (Jx) — a{Sx) is the induced moment 
(in units of gfis)- The single- ion mf Hamiltonian for 
each sublattice is then 



3 



A 
2" 



7 

-1 



(10) 



where 7'= 7(J), 7=^^- The mf- energies and states 
are then given by 



e± -±-[l + 7'= 

wis) 



|+; =u\e) 
\-)=u\e) 



(11) 



Here u=cos0, v=sin0 determine the rotation to the 
|±) mf eigenstates with 29 = tan^^(7'), furthermore 
sin 26* = (1 + 7'^)^ and cos 26* = 7'(1 + 7'^)^ where the 
convention 9 — Oa — —9b has been used. From the mf- 
single ion partition function Z=cosh/3e one obtains the 
mf- equations for the singlet-singlet splitting A'(T) in the 
ordered state 



A(r) 
A(r) 



A'(r) 



1 



(^4 



:etanh[(^)A(T)] 



)-(l- 
2^{S) 



(12) 



Here (Sz) = (S) is the pseudospin expectation value 
which can be replaced by ^ since A, furthermore 
A(0) — ^aj — ^. The control parameter of the singlet- 
singlet system is 



Therefore when ^ > 1 or T^v ^ A the saturation mo- 
ment becomes exponentially small. This distinguishes 
an induced moment system from a conventional magnet 
with degenerate ground state where the saturation mo- 
ment is a constant independent of Tjv. Note that the 
above relations also hold in the ferromagnetic case if we 
replace Je ^J//(0)= JAyi(0)+JAB(0). The temperature 
dependence of the induced moment is shown in Fig. ^ 
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(13) 



which characterizes the strength of exchange vs. size 
of CEF splitting. It determines the AF transition tem- 
perature where A(r;v)=0 as 



A 



2tanh~^(|) 2tanh~^(^) 
A 



— = tanh ■ „ 



(14) 



The staggered induced moment as a function of tem- 
perature is given by 



{J)t 



{Jh 



1) 



{J)o = a{S^h = \o,^{e 



4A^(T)-l]t 
-l)i 



(15) 



where {J)o is the saturation moment. To obtain a 
spontaneous induced moment (J)o > the control pa- 
rameter ^ must fulfil ^ > $c = 1- For ^ H- (5 ((5 < 1) 
slightly above its critical value = 1 the saturation mo- 
ment is given by 



2exp(-A) 



(16) 



FIG. 2. Temperature dependence of magnetic exciton soft 
mode frequency tj_(0)T (thin line) and normalized induced 
AF moment {J)t/{J)o (thick line) as given by Eq.(|l^ The 
small finite u!-{0)t below Tjv is due to exchange anisotropy. 
Physical parameters are the same as in Fig. 4. 



B. Magnetic exciton dispersion and AF soft mode 

The local CEF excitation at site i with energy A can 
propagate through the effect of the inter- site exchange 
whose transverse parts lead to transitions |e)i \g)i and 



l5>. 



simultaneously. This causes a widening of A 



into a dispersive band of 'magnetic excitons'. It should 
be stressed that these elementary excitations do not re- 
quire the presence of a magnetic order parameter as in 
the case of spin waves. Magnetic excitons exist already in 
the paramagnetic (PM) phase. As function of tempera- 
ture they may exhibit softening at a certain wave vector 
which signifies a magnetic transition that leads to the 
appearance of induccd_moments. This mechanism is well 
studied for Pr metalo and some compounds where the 
ordering wave vector may be incommensurate in contrast 
to the commensurate AF wave vector Q in UPd2Al3 . 
In the following a detailed description of the magnetic 
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exciton modes for UPd2Al3 is given since they are the 
bosonic excitations thought to mediate superconductiv- 
ity in this compound (Sect. IV. B). The calculation is 
done most conveniently in the linear response formalism 
in RPA which starts from-the dynamical susceptibility of 
the singlet-singlet systerr£3. 



(17) 



where i,j= lattice site, X, ii= AF sublattices (A,B) and 
a,/3= x,y,z are cartesian coordinates. Its Fourier trans- 
form is given in RPA by 



xiq,oj) = [l-u{uj)J{cf)] ^m(cj) 



(18) 



Here x, and J are 4x4 tensors, e.g. X—{x}\a.fif3 
with A, fj,= A,B and a, /3=x,y transverse cartesian coordi- 
nates. The tensor u is diagonal in A, fi and J is diagonal 
in a,/3. The expression for u(u;) is given in App.(B). To 
keep it in a simple form in the ordered state ((J) > 0) a 
real space back rotation by an angle 6^= d^= -0^ has to 
be applied to compensate for the effect of the molecular 
field. Under this rotation the exchange tensors transform 
as 



j'iXfi) = D{9^)J{Xfi)D^{e';:) 



(19) 



The magnetic exciton modes are then obtained as poles 
of the dynamical susceptibility in Eq.(p^) leading to the 
secular equation det x(?, w)=0 or 



\u'-\i,)-j'{q)\=Q 



(20) 



With u' ^{lu) given in Eq.(B5). The magnetic exciton 



branches obtained from the secular equation above will 
now be discussed for several cases. For brevity the primes 
in the transformed tensors will be suppressed unless ex- 
plicitly needed. 



1. paramagnetic phase, isotropic exchange 

In this case there are no sublattices A,B and only one 
exchange function J{q) — J oiq) + J n{q) appears, fur- 
thermore A' = A holds and Eq.(|ol) then yields the PM 
exciton dispersion 



uj{q) 



a^S)J{q) 



..(,-) = A[l-^^tanh^A] 



(21) 



At the AF point Q (PM zone boundary) we have 
J(Q)= Jd{Q)+Jn{Q)= Je and therefore 



w(g,T) = A[l-Ctanh^A] 



(22) 



This zone boundary mode in the paramagnetic BZ 
(PMBZ) becomes soft when T^ Tjv from above accord- 
ing to Eq.(p|). 



(23) 



^(4r) = i(A)2(^_ 
z In 4 



It approaches zero linearly above T^r. Below T^r the 
induced moment of Eq. ( |l5| ) appears staggered along the 
c- axis leading to two AF sublattices A,B. Then Q = 
(0, 0, tt) becomes the new center of the antifcrromagnetic 
BZ (AFBZ). One has the connection q' = q — Q, where 
q e PMBZ {\q^\ < tt) and ffj G AFBZ (|g^| < f ) as 
illustrated in the inset of Fig. H. The prime for the AFBZ 
will be suppressed i.f. unless explicitly necessary. 



2. AF induced moment phase, isotropic exchange 

For the solution of Eq. ( pO| ) one now needs the rotated 
exchange tensor in Eq. (p^) . For isotropic J(A/i) one ob- 
tains 



J'ix^l) = j{x^i)D{e^)D^{e'^) 
= J Dig) 



(24) 



J'Niq) 



1 



:JNiq) 



Where J'^iq), J'^i'l) intra- and inter- sublat- 

tice exchange tensors in the rotated (x',y,z') coordinate 
system respectively. From Eq. (po|) one then obtains the 
magnetic exciton dispersion 



[A' ~a^S){Jn{q)TrJNm 
[A'~a^{S){JDiq)TJNim 

1 - 7'2 2-e 



(25) 



1 + 7' 



72 



Here the last identity holds because T^ A where 7'^= 
A^ — 1 ~ — 1. Equivalently this may be written 

4(g1 = [A' - a'{S)J±{q)r [a^ {S)fi' Jl{q) 

JN{q) - (1 + i''r^JN{q); J±{q) - Joiq) ± JN{q) (26) 

In the AFBZ Lo^{q) and u>j^{q) correspond to acous- 
tic (A) and optical modes (O) respectively. Due to the 
property Jjv('f ± Q)=-JAr(9) the optic mode ^+(9) with 
q E AFBZ is the downfolded (shifted by an AF reciprocal 
lattice vector Q) acoustic mode ^^~{q) with q S PMBZ. 
This connection is illustrated in Fig. g (left panel). It is 
instructive to consider two limiting cases for this disper- 
sion: 



(i) e > 1 (r-^ 1) 

In this marginally critical case r>l and Eq 
to 



5h reduces 
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u:l{q) = A' - a^S)JD{q) T a'{S)JN{q) 



(27) 



The acoustic mode LO^{q) evolves from the paramag- 
netic soft mode uj{q) of Eq.(pl|) and the optic mode LO+{q) 
is obtained by downfolding lj- (q) . Therefore we have two 
strongly split A,0 modes in the AFBZ in the present case. 
Including only n.n. exchange along the c-axis in Eq.(|^) 
leads to J_d=0 and 3pf{qj—23ocosqz and we obtain with 
a^{S) =£_A~A: 

LO^{q) = A{1± cos q,) (28) 

The two modes are periodic in the reduced AFBZ. 

(ii) ^ 1 (r-^ -1) In this limit the magnetic excitations 
are exchange dominated and the influence of the CEF is 
negligible, therefore they can essentially be viewed as AF 
spin waves. One obtains from Eq.(25): 



«2(5))[4(0)-4(q)]1 



(29) 



where we used A' = = a^{S)oJe = a^{S)JNiO). 
Again, for only n.n. exchange along the c- axis {Jd—0, 
JAr(^=zJo cos(72, z—2) one obtains 



WAF iq) ^ ^^{q) = a^{S)zJa \ sin q^ \ 



sing^l (30) 



which is the spinwave dispersion of a simple AB- 
antiferrromagnet. Due to isotropic exchange it is twofold 
degenerate by which it is clearly distinguished from the 
strongly split A-0 exciton modes in the marginally crit- 
ical (C > 1) case. An illustrative comparison of the two 
cases (^ > 1, ^ ;» 1) is shown in Fig. ^ where Aaf= 
a'^{S)zJo is the spin wave band width with AafI^ — 
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FIG. 3. left panel: Magnetic exciton branches u}± in the 
AFBZ (FAA' is equivalent to AFA, see also inset of Fig. 8. 
Physical parameters as in Fig. 4. right panel: Simplified model 
with only n.n exchange along c. Comparison of magnetic 
exciton branches lu± for the CEF- dominated case (C >1) and 
(twofold degenerate) AF spin wave dispersions loaf for the 
exchange dominated case (^ S>1), normalization to A and 
Aaf of Eq.(bd) respectively 



3. AF phase, anisotropic exchange 

In the isotropic case a gapless Goldstone mode for q — > 
(AFBZ) will be present for any ^ >1 above the criti- 
cal value. However due to the uniaxial symmetry of the 
effective exchange tensor Eq.(^ in principle one has to 
expect gapped magnetic excitations. For simplicity we 
assume that the anisotropy is independent of the spin 
pairs, then we may define uniaxial Jf) N{q)i Jh N^'t) 
change functions in obvious notation where a,c refers to 
the hexagonal directions. The magnetic exciton modes 
may again be obtained as solutions of Eq. (po|) , however 

the transformed J'jj [q) and J'^ {q) are no longer propor- 
tional to the unit matrix as in Eq. (p^) . Instead one has 
for the restricted (x',y) coordinate system: 



J'd = D{ 



9^) 



ulJf^+vlJ% 
J% 



„2 ja _ 2 T. 



N 







ja 



(31) 



(32) 



where u^ = cos 9r, Vr = sin 6r and the rotation angles 
for the two sublattices are Of = 9r = —29 and 6^ = 
-Or = 20. Therefore (App. B) u^ = cos26' = a = 
(1 +7'^)"^ and v^ = sin26' = 6 = 7'a. From Eq.(|20l) and 
using Eq. ( |3l| , |3^ ) we then obtain the magnetic exciton 
modes in the most general case considered here as 

Loliq) = [A' - a^{S)Mq)Y {a^S)f^'^ J{q)l 

J±{q) = [Jd ± J%1 + (1 + i^)-\Jf,' ± Jn) 

J±{q)^{l + i^r\jfi^±JN) (33) 
1 



Jd. 



N 



(Jd. 



N 



J- 



'D,N 



D.N 



D,N) 



'D.N) 



In the case J^= J^^ and 3% this reduces to the 



isotropic result of Eqs.(^23,2' 



C. Model calculation, comparison with INS results 

The simple AF structure was determined in Ref. [ |2^] 
and subsequently inelastic neutron scattering (INS) ex^ 
periments on single crystalline UPd2Al3 were performedLl 
in an energy range up to 10 meV in the whole AFBZ 
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to determine the collective magnetic excitations. It was 
concluded that for q\\ c*-axis well defined magnetic exci- 
tation modes exist. For q _L c*-axis their dispersion was 
also determined, however line widths are generally larger 
in this case. An extrinsic damping mechanism via cou- 
pling to itinerant 5f- electrons was concluded from the 
absence of a zone center gaxL pHowever later more accu- 
rate low energy experimentsEj'tiJ have shown that a gap of 
1 meV exists at the AF wave vector (AFBZ zone center). 
Within the present RPA theory the exchange anisotropy 
introduced before can only generate small gaps < 0.2 
meV because of the vicinity to the quantum critical point 
at ^c=l of the AF singlet- singlet system as discussed be- 
fore. The observed zone center gap may have a dynami- 
cal origin due to terms neglected in RPA. In addition the 
q ~0 low energy excitation modes will be strongly shifted 
and broadened due to the interaction with 5f- conduction 
electrons as proposed and investigated in Ref. ( |l^) . 




r M K r A A' 



FIG. 4. Magnetic exciton dispersion of UPd2Al3 in the 
hexagonal AFBZ. Squares are data from Ref. [9]. Solid line: 
Calculation from Eq.(p^ using the exchange parameters Ji 
(i=0-4) in meV Jo= -1.63, Ji= -0.4, Ja^ -0.31, Ja^ 1.75 and 
J4= 0.6 as defined in Fig. 1. For technical reasons an a-c 
exchange anisotropy of 1.1 is assumed which leads to a very 
small F- point gap. 

In the present discussion the focus is on the funda- 
mental origin of the superconducting pair potential and 
the resulting gap function within the weak coupling ap- 
proach. For this purpose the details of the low energy 
spectral behaviour of the soft mode (AFBZ) is not of 
primary interest. It is more important to obtain a good 
global description of collective magnetic excitations for 



all wave vectors and energies since it is the total spec- 
trum which determines the strength and anisotropy of the 
pairing potential. For this reason we use the RPA theory 
developed in the previous section to describe the mag- 
netic excitations found in Ref.( ||) for the whole AFBZ. 
Going beyond RPA for the magnetic exciton spectrum 
would also necessitate a strong coupling approach to the 
pairing problem to stay on the same level of theoretical 
treatment for magnetic and superconducting sub- sys- 
tems. This challenging problem will be left for a future 
investigation. 

We will now discuss the experimental and theoretical 
dispersion of magnetic excitons for UPd2Al3 unfolded 
along typical symmetry lines of the hexagonal AFBZ. 
For the singlet- singlet system the magnetic exciton dis- 
persion according to Eq. ( ^3| ) requires the following quan- 
tities: CEF- gap A, exchange constants for the various 
neighbor shells and the uniaxial anisotropy constant. For 
the CEF splitting we take A= 6 meV which corresponds 
roughly to the center of the magnetic exciton band as 
seen e.g. in the simple model of Sect. (III. B. 2). Further- 
more this is also approximately the value obtained frMSji 
fitting the susceptibility as function of temperaturell3. 
Since Tn (14.3 K) < A (66 K) we are in the CEF dom- 
inated regime: According to Eq.(p^ ^ — (tanh 2^)^^ — 
1.015 is only slightly above the value for the QCP at 
^c=l which leads to a saturation moment {J)o ^ 
and hence A' ~ A. Therefore the excitation spectrum 
will be almost like the magnetic excitons in the param- 
agnetic regime. The exchange model used for the local 
moments of UPd2Al3 is illustrated in Fig. 0. A consider- 
able number of neighbor shells has to be included to get 
an agreeable fit to the experimental data. The result is 
shown in Fig. || with the exchange constants in the figure 
caption corresponding to the notation in Fig. ^ The cal- 
culated and experimental dispersion are strongest along 
the hexagonal c-axis where the modes are also well de- 
fined i.e. the line width is small. In addition it shows a 
pronounced dip at the K-point. Because of the relation 
LOj^{q ± Q) = Lo^{q) only uJ-{q) is plotted. In the sector 
FA it corresponds to the acoustic mode and in the sector 
AA' it is the optical mode, this is illustrated in Fig. ^ 
(left panel). The large mode energy at A' gives direct 
evidence that one has a CEF- dominated singlet- singlet 
system; in the opposite exchange dominated regime with 
hypothetical T^r ^ A one would have almost degenerate 
acoustic and optical branches close to u>af shown in Fig. 
^ (right panel) and in this case the excitation energy at 
A' should be very small again. It is not entirely clear 
how realistic the dip at the K- point reported in Ref. ( |^) 
is in view of the fact that the line width becomes much 
bigger than the mode energy at this symmetry point. 
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IV. THE SUPERCONDUCTING PAIR 
POTENTIAL 

In the known superconductors the formation of Cooper 
pairs is thought to be mediated by the exchange of 
bosonic excitations of the medium. This is a simph- 
fied picture of a comphcated many body problem. More 
generally the pair formation is the result of a singular 
behaviour in the low energy part of the two particle scat- 
tering matrix. In the simplified picture this singularity 
is described by the repeated electron-electron scattering 
via boson exchange. In a nonretarded weak coupling ap- 
proach this is equivalent to the formation of a bound 
state in the appropriate scattering channel. 

Two origins of Cooper pair formation are commonly in- 
volved for describing real superconductors: (i) the (BCS) 
electron-phonon mechanism in conventional supercon- 
ductors where the effective attraction needed for bound 
state formation is provided by the exchange of a real 
frequency (propagating mode) phonon. It is spin inde- 
pendent and allows only for the formation of s- wave spin 
singlet Cooper pairs, (ii) the spin fluctuation mechanism 
in strongly correlated electron systems (HF-compounds, 
cuprates etc.) where the Cooper pair formation is due to 
exchange of an imaginary frequency spin fluctuation of 
conduction electrons, predominantly close to an AF zone 
boundary wave vector (overdamped " anti-paramagnon" 
mode). The latter mechanism which is in principle rel- 
evant also for UPd2Al3 , has originally been intoduced 
in the context of superfluid He^ by Nakajiman3 in its 
paramagnon form. Later it has been adopted for vari- 
ous Heavy Fermion supercoductors by many authors, see 
e.g. Ref.( 23) because it naturally leads to the possibility 
of unconventional pair states thought to occur in these 
compounds. It involves the spin degrees of freedom in 
a rotationally invariant way and allows for spin singlet 
and triplet pair formation, (iiil As a novel mechanism 
it has recently been proposedEZI that in U- compounds 
like UPd2Al3 with 5f- electrons of dual character the ex- 
change of magnetic excitons is at the origin of supercon- 
ductivity. This mechanism is distinctly different from 
those discussed before. On one hand it is mediated by a 
real frequency propagating mode, the magnetic excitons, 
on the other hand, as in the previous case it involves the 
spin variables, but due to the presence of a CEF splitting 
the rotational symmetry is broken. 

In this part of the present work we will derive and dis- 
cuss the nonretarded effective pair potential of this new 
mechanism in detail. Because it is important to under- 
stand its differences as compared to the spin fluctuation 
mechanism usually invoked for HF metals we flrst give a 
short summary of the salient ingredients of the latter. 



A. The spin fluctuation model of pairing 

This pairing mechanism may be dominant in many 
strongly correlated electron systems where one has a 
screened on-site Coulomb interaction of heavy quasipar- 
ticles described by 



(firni (r)r7,| (r) 
V ^ feT -fe 



y ^ k\ -fei 
fe,fe ' 



(34) 



Sumifting appropriate two particle scattering diagrams 
in RPAl3 Eq.(|3j) leads to an effective pair Hamiltonian 
{q = k' -k) 



kk' ,Q.(3^5 
ka —k'y —k'S k' (3 



(35) 



where the sum in parentheses describes the effective 
two particle interaction and tr, a ' denote Pauli matrices. 
It contains a spin- independent potential scattering of 
strength p{q) and an exchange scattering of strength 
y s{<i) which are obtained in RPA as 



Vp{q) = 



It 

2^c 



VM = 



2 c 



1 - icXoiq) 



(36) 



Using the projectors for pair spin singlet (S=0) and 
triplet (S=l) channels respectively. 



^0 

Pi 



(37) 



with the Kronecker product {aa ')af3:-yS— f?a/3f? ^5, one 
can separate ^eff{q) into the irreducible potentials of 
singlet and triplet Cooper pairing: 



VM = VM - ^Vs{q) 0, -Ilxiq) 



VM^VM + VM 



(38) 



where the last approximation in the above equations 
holds for wave vectors with enhanced spin fluctuations, 
i.e. when x{<t) > Xo(<z)- Here 



xo{q) 



1 - icxM 



(39) 



is the RPA susceptibility of conduction electrons which 
has its maximum at a nesting vector Q which may be 
the commensurate AF vector. The effective potential 
may be written, using spin projectors Pq, Pi and defining 
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the sequence of states as (ff , t|., |.t, J,J,), in the following 
form: 



Veffiq) = Vo{q)Po + Vi{q)Pi = 














(40) 



where 



to Eq.M 



according 

Note that the total spin component Sz = ^{crz + cr'^) is 
conserved due to rotational invariance of Ve/ / , therefore 
matrix elements corresponding to a change of the z- com- 
ponent of the pair spin like (tt \Veff\ il) = (ii \Veff\ 11) 
vanish. One may explicitly check by diagonalisation that 
"Veffiq) has singlet |0) and triplet |1, Sz = 0, ±1) eigen- 
states with energies Vo and Vi respectively. The mo- 
mentum dependence of xio) finally determines which or- 
bital pair state will be realized as disussed in Ref. [ |2^] 
who conclude that even parity singlet pairs are favored 
in general. 



B. Pair potential due to magnetic exciton exchange 

In the dual model of Eq. (|^) the on-site Coulomb inter- 
action of the itinerant 5f- electrons is neglected and only 
indirect interaction via intermediate CEF excitations due 
to the last term is considered. Since the latter can prop- 
agate this means an effective conduction electron inter- 
action via the exchange of magnetic excitons. This will 
now be derived by similar diagrammatic methods as in 
Ref. [ , however the details will be quite different for 
the present case. In this section we neglect the main ef- 
fect of AF order on the pair potential and consider only 
the paramagnetic case. For the derivation of the pair po- 
tential the Hamiltonian Eq.(|l|) will be first rewritten in 
terms of more convenient bosonic variables for the CEF 
singlet-singlet excitations introduced via the pseudo spin 
Holstein- Primakoff representation 



St^ = a,- 5.- —at Sf = — — a. 



(41) 



Here local boson that describes CF- excitations 

\g) — > |e) at site i. For T^ A the influence of unphysical 
states (ai)"|0) with n>2 is negligible. This condition is 
very well fulfilled since T< < A. Then the on- site 
exchange 

H,f = -2/0(5 - 1) Et^--^- + 'tJr + s-Jt] (42) 



can be transformed, using — e(^ — S^) and 
aS"^ and Eq.dfl) to obtain 



fei k-k' k'l k'-kl 



kk' 



-^'S^(ci CTt- — ct cr, )a^' ~ apt (43) 

k'l- feT k'l kl' n + k-k' « ^ ' 



kk' q 

where I=a(g-l)Io and Iz=e(g-l)Io. For a, e see Eq.(||), 
furthermore 



iqRi 



(44) 



creates CEF- bosons of momentum q. The total Hamil- 
tonian of Eq. (|l|) is then given by 



k<T 



-k 1- 




-k' ^^ 



kH 



k1^ 



-k ^^ 





-k' if 



k"^ 



FIG. 5. Lowest order (~I*) polarisation diagram (above) 
with q = k' + k and exchange diagram (below) with q — k' — k 
contributing to the scattering potentials Va and Vd of Eq.(|5c|) 
respectively. The zig-zag line represents the boson propaga- 
tor Do and the bubble is equal to the conduction electron 
susceptibility Xoio) in the nonretarded approximation. Small 
arrows represent conduction electrons spin. Note that the up- 
per diagram describes a process which does not conserve the 
total pair spin component (s^-l-s'^) leading to a pair poten- 
tial Eq.(p2[) which is not rotationally invariant in spin space. 
Summation of diagrams of odd order in leads to Vj, and 

Ve. 



The interaction term in Eq. (^3|) contains both 1- boson 
absorption or creation parts as well as 2- boson scattering 
parts. The latter describes conduction electron scattering 
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from quantum fluctuations in the occupation of the two 
singlet CEF- states. Its contribution will be neglected in 
the following discussion. The boson propagator is given 
by 



%(r) = D+{-t) = -{T,a,iT)a+m 
and its Fourier transform (oj, 

D^{iuj„) = Di{-iuj„) 



q\ - ,-~q 

2n7rT) by 
1 



A 



(46) 



(47) 



The momentum q for the noninteracting CEF- bosons 
is a dummy variable. Note that in the zero frequency 
limit (D^Di)i^^^o = sT The lowest order nontrivial po- 
larisation and exchange diagrams for the effective e-e in- 
teraction vertex due to the 1- boson part of Eq. ( ^3| ) are of 
order ^I"* and are shown in Fig. |5[ For the nonretarded 
e-e vertex function Fq {q) that enters the effective pairing 
Hamiltonian we only need the zero frequency limit of Fig. 
H in which the diagrams lead to equal contributions: 



+ ^ ^ 14(17) [ct - c T/ 1 cr,-h + ct - c rit-CTn] (51) 

kk' 

+ Vrf((7)fct - C T/ I CT/^ + ct - C Tz-hCr/il 

av^n ^1 -k'l fc't fcj^ -fc't fc'xJ 

kk' 

Defining the effective e-e interaction potential by 
Heff = \ ^ E ^^pU^'^l.Ak.'-k'sCk'p (52) 

kk' ^afi^S 

we get the explicit effective interaction matrix by com- 
parison with the previous equation, using V±=Va±Vf, 
and again defining the state sequence as (tti Ti; it, ii)- 



1 

















v+) 
-v+ 









v+ 







(53) 



-Xo{q,0)? (48) 



Here the noninteracting susceptibhlity Xoio) = 
Xxx{<t) ~ Xyyio) is given by the Lindhard function 



Xoiq) 



fk f k+q 



(49) 



The summation of polarisation and exchange diagrams 
as in Fig. || to infinite order in RPA leads to the follow- 
ing nonretarded e-e interaction potential Yap-ysi't) for the 
various channels of initial {(3S) and final (aj) spin con- 
figurations: 



Va 










1 




i-n 2 ^A(i-ro) 


A(i + ro)^ 


Vb 




(50) 




I'Do 1 , 1 


1 




i-rl 2 ^A(i-ro)^ 


A(i + ro)^ 


Vc 













Here Va(g), Vb(g) with momentum transfer q = k + k' 
correspond to the polarisation and Vc((?), Ydiq) with q = 
k — k' to the exchange diagrams of Fig. ^. Adding up 
all contributions and transforming fc — )■ — fc in the first 
two (a,b) terms one obtains the effective e-e interaction 
Hamiltoninan where now q = k — k': 



^eff - 



kk 



Comparison with Eq. ( [40|) shows that interaction com- 
ponents V|||| etc. are absent due to the fact that bo- 
son absorption or emission always leads to spin flips. 
On the other hand there are additional components like 
=Va-l-Vb=V+ which do not conserve the total pair 



Y. 



TTii 



spin component s* ^(s^-l-s'^). Thus magnetic exciton ex- 
change does not lead to a rotationally invariant effective 
pairing potential as in the case of the spin fluctuation 
mechanism in Sect. (IV. A). This is a direct consequence 
of the presence of a CEF splitting. The decomposition of 
Yef / into irreducible contributions as previously given in 
Eqs.(^0|) is therefore more involved and will be discussed 
below. 

The static e-e vertex function rQ{q)= ^o{q)- 
J2q^oi'l) is related to the Fourier transform of the 
effective exchange interaction Jrkky (q) of pseudo spins 
via the the relations 



{S) 



jRKKviq) 



(54) 



where Xo(^) was given in Eq. (^,p9|) . In addition the 
pseudo spins have a superexchange iff{q) whose con- 
tribution in static RPA limit is described by equivalent 
diagrams thus leading to a total vertex function 

nq)^T',{q)+Tff{q)^ 



[J. 



RKKY 



[<f) + Jjj{<f)]^^J{^ (55) 



This result may now be substituted in Eq.(^0|). As 
mentioned in Sect. (II) the total exchange function J((j) 
will be parametrized by comparison with the exper- 
imental exciton dispersion. It has intra- and inter- 
sublattice parts Jd(<?) and JAf(q) respectively with 
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J(g)=JD((7)+Jjv('z)- Absorbing the constant jf^^^oio) 
into a renormalized CEF splitting A and replacing A 
A' for the AF ordered case we then obtain, defining 
ujRiq) ^ 2A' ~UJ^{q) 

A'[l - rm - A' - a^S)JDiq) - a^{S)JNiq) = ^-i<l) 

(56) 

A'[l + Tiq}] - A' + a^S)JD{q) + «'(^) Jjv(<?1 = LUR{q) 

Here uj-iq) is nothing but the excitonic dispersion 
branch which becomes soft at the zone boundary Q of 
the PMBZ at Tjv and evolves into the acoustic branch 
with its minimum at the new T- point of the AFBZ as ex- 
plained in Sect. (Ill) Note that the energy ujR{q)=2N — 
uo- (q) does not correspond to an excitation branch. Since 
2A' ^ t^-io) for all realistic cases u!R{q) will never be 
small and a;_R(g)"^ is regular throughout the AFBZ. Us- 
ing Eq.(p6|) the effective potentials of Eq.(58) can now be 
written as 



V„=--r 



1 



1 



2 LO^{q) LUniq) 



1 



2 LU^{cf) ujRiq) 



(57) 



Inserting this result into the interaction matrix of 
Eq.(p3|) the possible irreducible pair potentials (i.e. the 
eigenvalues V^, k=0,u,v,w of Vejf) are given by 



Vu = -{Va+Vb)=P 



1 



w_(g) 



Vy = (Va + Vb) 
K„ = {Va + Vb) 



-P 
-P 



1 



1 



(58) 



with the correspondence Vo=V_ , Vu= -Vv= -V^= V+ 
to Eq . (|53|) . The pair states obtained as eigenstates from 
Eq.(p3|) are classified according to their Sz=ti (=±5) 
quantum numbers. Later on when the effect of AF order 
is considered, a representation with Sj;=-|-,- (=±i) may 
be useful. It is obtained by |-|-) = Ry{^)\ T) etc. where 
Ry = -i=(l — iay) is a rotation by ^ around the y- axis. 
Therefore we also give the relation between the two no- 
tations for the eigenstates ipK, or ip^ which correspond to 
Vk, K=0,u,v,w: 



V'o = ^0 = ^(1 Ti) - 



it)) 

= -C = ^(1 TT) - I ID) 
ID) 
it)) 



i^v = C = ^(1 TT) + 



(59) 



= C - ^(1 Ti) + 



Here ipo describes singlet pairs. In the spin fluctua- 
tion model (■0K,K= u,v,w) correspond to the degenerate 
triplet pairs. In the present magnetic exciton model ac- 
cording to Eq.(|5^) obviously this degeneracy is partly 
lifted because the basic Hamiltonian Eq. ( ^ , ^5| ) is not 
rotationally invariant in spin space as a result of the cou- 
pling to CEF split localized states. The triplet is split 
into a nondegenerate state ipu and a doublet (ipvj^'w)- 
As a consequence in the present case the irreducible pro- 
jector representation for the pair potential like Eq.(p3|) is 
different. The projectors to the eigenvectors which sat- 
isfy PiPj =Pi5ij may be written as 

Po = PS = ^{l-aa') 

Pv = p: = \{i + ^.4 - + <yrA) (60) 

Pn, = Pu = \{l + <yx(T', + OyO'y - <7,a'J 

The projector to the doublet subspace P d=(Pv+Pw) 
and to the total triplet subspace Pi = (P,u+P^,+Ptu) are 
given by 

Pd = 1{1 + ct..K) 
Pi = 1(3 + <?<?') 



(61) 



Then similar as in Eq.(pO|) Ve//(g) as mediated by 
magnetic excitons (Eq. (p3[) )may finally be written as 



Veffiq) 



:Po 



Pu 



ujR{q) uj-{q) uj-{q) 



-{Pv+Pw) 



(62) 



There is an intuitively simple interpretation for the 
CEF- induced triplet splitting present in this model: The 
conduction electrons can scatter from virtual CEF- exci- 
tations A only with their s^; (or Sj,)- component, therefore 
the pair states will experience a potential ~ (/^/A)(s^°*)^ 
in lowest order which splits the (jp% , k=u,v,w) triplet into 
nondegenerate -ipw (Sx°*=0) and -tpu '4'v = ^l) de- 
generate pair states. 

The previously developed theory of magnetic excitons 
which explains the experimental observations as shown 
in Fig. ^ can now be used to calculate the pair potential 
from Eqs.( p8| , |62[ ). However only the q^;, qj^- averaged 
potential is needed for the gap equations as explained in 
the next section 



V. THE SUPERCONDUCTING GAP 
EQUATIONS 

The effective e-e interaction may lead to a supercon- 
ducting state with an order parameter given by the gap 
matrix (q = k — k') 
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(63) 



which satisfies the m.f. gap equation 

- A*"'- 1 
Aj- = - V V,''Z,t tanh 

(a/3)(7'5) 2K- 2 



Here we defined V^^^^j'j^^^) = V^J^ and (a/3), (7-5) 
TT,Ti,iT,ii arc double spin indices in the notation for 
s^- quantisation and E^: are the superconducting quasi- 
particle energies. The gap matrix may be decomposed in 
the canonical way as 



(64) 



A* 



3 



(65) 



where d§, d^=(d^,df ,d£) correspond to the pair states 
■00, (■i/'u,V'ti, "01)) respectively introduced before. Explicitly 
we have the relations in both s^- and s^, quantisation 
notation: 

4 = i(Af, -AfT) = i(Ai_-A^:+) 



< = -2(Ar 



Ati) = 2(A 



At,) = ^(Al 
AfT) = -^(A* 



+ A'::_ 



(66) 



-) 



Using the projector representation Eq.(^2|) for the ef- 
fective potential and the identity 



(67) 



where A'^={A(ct^)} is a four component vector acted 
on by the 4x4 projector matrices of Eq.(|6^). We then 
obtain the scalar irreducible gap equations (k=0,u,v,w) 

di = -Y^vy^'F^^,di 





4- 


^ COS kz ■ 


^i.(rf) 




= (4,0,0)^ 


■• X sin kz 




d-k 


= (0,4,0)- 


^ y sin kz 


Bsu i^i ) 


4 


= (0,0,4)- 


■> z sin kz 


Ain(rr) 



(69) 



These representations have even (dg) and odd {dj^, k — 
u^v,w) parity respectively The basis functions have all 
odd character with respect to the transformation k 
k±Q, namely 4''''^= "4 ('^= 0,u,v,w). In 1)211 there is 
no symmetry reason why (dj,dj^) should be degenerate. 
In fact we show later that the AF order connected with 
symmetry will lead to a splitting into nondegenerate 
d^ and d^J. superconducting states. 



Solution of gap equations for cylindrical 
symmetry 



The explicit solution of Eq.(|6q) using the potentials 
of Eq.(|5^) is greatly facilitated by taking into account 
the corrugated cylindrical FS shape modeled by Eq.(||) 
and shown schematically in the inset of Fig. ||. In this 
case Eq.(^8|) can be approximated by a one- dimensional 
integral equation in the variable kz- Explicitly then 



dK{kz 



E 



Vi{kz - k'^)d^{k'^] 



tanh|£;(fc^,/c^) 



E{k'^,k'z 



(70) 



where k±={'kx,^y), kx=^±cos(l), 'ky=k±sm(j) and ^ is 
the azimuthal angle in the ab- plane. For the above sep- 
aration of variables it is assumed that the 0- dependence 
of the order parameter can be neglected, this means we 
restrict to representations of the type given in Eq.(|69|) 
and higher harmonics. In addition for the cylindrical FS 
sheet in UPd2Al3 e{k±,a) depends only on the modulus 
kj_. Therefore only the q2:,qj,- averaged pair potentials 
Vk,{<Iz) appear in the gap equation where we define 

^i^ilz) = -ri I dqj,dqyVi^{qx,qy,qz) (71) 



Ft 



k 2 ^ 2 ^ 



(68) 



Ek = [i^k-^^f 



Before solving these equations explicitly we note that 
the gap functions can be classified according to the crys- 
tal symmetry group which in this case is the orthorhom- 
bic group D2/1 due to the presence of the AF order pa- 
rameter with ffi^ II X. The corresponding irreducible 
representations of gap functions can be listed as 



Using the expressions in Eq.(58) and the calculated ex- 
citon dispersion for UPd2Al3 in Sect. (Ill) we obtain these 
averaged pair potentials as shown in Fig. ^. Because of 
the singular behaviour of {q) at the AF- point Q the 
absolute value of averaged potentials |Vk;((7z)| (k=u,v,w) 
show a pronounced maximum at (5=(0,0,7r). On the 
other hand the averaged singlet pair potential Vb(q2) is 
rather smooth with a fiat maximum for qz around 0.47r. 
This is due to the nonsingular nature of Yo{q) as dis- 
cussed in Sect. (IV). 
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FIG. 6. qz- dependence of the absolute valued of q^,, q^- 
averaged pairing potentials Vk in the PMBZ for even (Vo > 0) 
and odd {Vu > and Vv, Vm < 0) parity of the gap function. 
qz=0 corresponds to the A'- point and qz = f- to the F point 
of the AFBZ in Fig. 4 respectively. 

The gap equation Eq. ( [70| ) may now be completely fac- 
torized by expanding VkJ^z) (qz=k^-k'2) in lattice har- 
monics according to 

m 

Vi^{cosmkz cosmk'^ — s'mmkz sinmfc^) (72) 

m 

Furthermore we may expand the even (k— 0) and odd 
parity {k— u,v,w) gap functions into 

do{k^) = ^ d^cosmfc^; df,{kz) = ^ d^sinmfc^ (73) 

m>0 m>l 

The order parameter may now be in general repre- 
sented as vectors do — {d™} ^^"^ = {d™} and the 
gap equations for the components obtained from insert- 
ing Eqs.(|7|,^) in Eq.(|70|) finally can be written as 



l = a,iV(0)FrG™K 
1 



aUK)F[T,K,d^]dk'^ (74) 



r - . r- tanhi/3S(C,fc:) 

Here w^L is the magnetic exciton band cutoff energy, 
CTk=-1 for K = and crK=l for = u,v,w. Note that ct^ 



is due to the different signs for even and odd parity parts 
in the harmonic decomposition in Eq.(|7^). Furthermore 
we have for k=0: 



i^(^,A:.) = [e' + (^Ccosnfc,)^]^ 

n 

amikz) = cos^{mkz) 
and for k,= u,v,w similar equations 

E{tkz) = [e + iY.d:smnkz)']i 

n 

a.m{kz) = sin^(mfc^) 

At the critical temperature Tc d^—Q and Eq.(^ 
plifies to 



(75) 



(76) 



Sim- 



F{T,) = Ar(0) 



G^[0\^a^F{T,) 



tanh \l3ci 
I 



d^ = iV(0) ln(A/3ct^l) 

(77) 



Here A=1.13 and is the average of am{kz) in the 
interval [0,7r]. One has am=l for k=0, m=0 and a™ ~ ^ 
for K= 0,u,v,w and m>l. This leads to the weak coupling 
BCS- formula for the transition temperature: 



1.13(^1) exp(- 



1 



a„7V(0)|y- 



r) 



(78) 



For T™ >0 the conditions Fq™ <0 for the even d™ (k = 
0, m > 0) singlet and V"^' >0 (k = u,v,w,m > 1) for odd 
parity d™- states must be satisfied. This difference in the 
required sign of the Fourier components can again be 
traced back to the decomposition in Eq.([72|). The Fourier 
components in Eq. ( [72[ ) with m= 0-4 are shown in Fig. |^. 
Within the weak coupling approach the sc phase realised 
is the one with the highest T™ , i.e. the highest >0 
(k=u,v,w) or —V(^ >0 {k=0) value. From the calculated 
values in Fig. ^ we notice that the degenerate V^T^^^ is 
the most favorable which is slightly larger than V"^"^^. 
The former corresponds to an odd parity doublet 



d{kz) = {d]j sink z,d]^ sink z) 



(79) 



with \d]j\ = \d]^\. It has node lines of the gap at k2=0. 
As mentioned before this degeneracy is accidental in the 
present model due to the neglect of AF ordering and 
therefore there is no symmetry determined phase relation 
beween the d\ and d^ amplitudes. The second slightly 
less favorable odd parity state d^ sin 2kz is nondegenerate 
and has line nodes at k^^O, ^ 

An important conclusion from this analysis is that for 
the present pure magnetic exciton model the singlet state 



do{kz) 



(ig COS kz 



(80) 



is not stable because Vj >0 is repulsive instead of at- 
tractive as required in this channel. This is due to the 
fact that the potential Vo{qz) in Fig. 0(dashed line) has 
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its maximum not at the PM zone boundary (q^ = tt) but 
rather close to qz = OAtt. This behaviour can be traced 



back to the form of Yo{<t) in Eq.(58) which shows that 
the denominator 2A' — {q) and not the exciton mode 
energy uj- (q) which leads to the maximum at q^ = tt for 
the odd parity state potentials Vk(9z) at qz = tt in the 
PMBZ. This result is robust against the details of the 
exciton exchange model. 



0.3 



> 
E 




-0.3 



m 



FIG. 7. Fourier components of Vk{(Iz) according to 
Eq.(|7^). dashed line (diamonds): k—Q, full lines: k=v,w 
(degenerate,triangles) and K=u(squares). The Fourier com- 
ponents with m<2 which lead to a nonzero TJ!^ in Eq.(|7^) 
have been designated by arrows. The odd parity (v,w) com- 
ponent (m=l) is the most favorable. 



B. Effect of AF order on superconducting gap 
equations 

In the previous treatment the influence of AF order 
on the pair potential was largely neglected. Sofar it only 
enters through the renormalized singlet-singlet splitting 
A' in the magnetic exciton dispersion which is of minor 
importance. A much more pronounced effect is due to 
the competition of sc pair formation and reconstruction 
of quasiparticle states at the PS by AP long range order. 
This leads to an important modification of the effective 
pair potential for wave vectors close to the AF ordering 
vector Q. This is also known from AP phonaiumedi- 
ated superconductors such as the borocarbidesEJS. The 
application of this idea to the present case is formally 
similar but more involved because of the more general 
possible pair states as compared to conventional electron- 



phonon singlet superconductors treated in the above ref- 
erences. In this section we show how AP order modifies 
the magnetic exciton pair potential of Eq.(|62[) and the 
effect on the sc states as discussed above. 



1. Transformation to magnetic Block states 



The AP order parameter 



(Ji) = (J)QXCOs{QRi) 



(81) 



leads to an additional periodic potential for the con- 
duction electrons due to He/ in Eq.(^) which has twice 
the period (2c) of the lattice potential (c). It is given by 



cl . cr.) (82) 



and leads to a reconstruction of the conduction band 
states (see e.g. Refs.[ ^,^) which may be described by 
a Bogoliubov transformation to the new magnetic Bloch 
states of the AP lattice given by 



at = a.-ct 
fci fcT 



■^fe k+Ql 



at ^ = B r ^ct 

k+Q2 ' -k-Q fc| 



-fc-Q k+Ql 



(83) 



The Ct^rr-,r 



ir are the real 

ik 



magnetic Bogoliubov 
coefficient^^lSa below T^r (not to be confused with Uj^jVj^ 
which transform to the superconducting quasiparticles 
below Tc). The effect of AP on the pair potential Eq.(|6^) 
is included by expressing the latter with the transformed 
pair operators Bj:={aj:^a j:^, , {aa') — 11.12,21,22} in- 
stead of the pair operators b^—{c^^c_^^,, {era') =tT;Ti 
, It, 11} in the paramagnetic state. This transformation 
b^=Tj:B^: is furnished via the symmetric 4x4- matrix 



( -I 

^kl^k 
^kPk 



-I 
"I 



Pi 
p\ 



'^kh 



«E ^kh P\ 



(84) 



Then the effective pairing Hamiltonian transforms into 
(q = fc'- fc) 



kk' 



\y^!Mh-\Y.^''^J^sf^^h (85) 



kk' 



using the projector representation of Eq. (|6^) the effec- 
tive pair potential may be written as 



(86) 
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As long as k, k' are far away from the AF Bragg planes 
k = ±5(5 one has T^:,Tg, ~ 1 and P„ ~ P„ stih projects 
to the eigenstates of Eq. (B9h . However close to the Bragg 



planes Pk 7^ Pk and Vef / (q) and hence the gap equations 
will be strongly modified. ExpHcitly we have 



Pv Pv ~l~ bj^./b^Ptu ~\~ Pvw H~ bj^Pujy 



(87) 



Here we defined the transfer operators Pvw— PLu— 
\'>pv){'ipw\ between different pair states, this means that 
AF order mixes ipy, 'ipw states for wave vectors close to 
the Bragg planes, contrary to -00 ' V'« which stay in the 
same subspace. The coefficients in Eq.(|87|) are related to 
the magnetic Bogoliubov coefficients via 



g/^g- i,_T;? M -pi) (88) 

Where l(^,)=lo(g-l)(J)o 

C. Modified gap equations in the AF state 

The above relations enable us to derive the modified sc 
gap equations in the presence of AF order. The effective 
pairing of Eqs. (^ , ^ ) leads to gap equations analogous 
to Eq.(|63|) for electrons in magnetic Bloch states: 



(89) 



Multiplikation with a projektor from the left and 
using the identities P«P«u,=P«u,, P^Pwi, =0 etc. we ar- 
rive at the modified gap equations for k=0,u: 



k' 



(90) 



This shows that for the V'o ^ i'u P^i^ states the effective 
pair potential in the presence of AF is given by 

V,{k ~ k') = (4 - /32)K(fc - fc')(«| - (91) 

which vanishes at the Bragg planes k = ±^Q, i.e. 
states which are connected by the AF ordering vector 
do not contribute to sc pairing. This leads to the well 



known depression but usually not destruction of super- 
conductivity by AF order, see e.g. Ref . ( |2^) . Far from the 
magnetic Bragg planes Vk(A; — k') ~ V^ik— k'). Further- 
more the gap function entering the quasiparticle energy 
E, is d^ which has an additional zero for k = i^Q due 
to the prefactor (a| — /3|). This factor does not change 

sign at k = i^Q and therefore does not lead to an ad- 
ditional node line. The modified gap equations can now 
be derived in a similar way for k=u,v: 



E 

k' 



ibj:Vy{q) 
VUq) 



d^' 
d^' 



Ft, 



(92) 



Obviously the two states (d^,d^) are mixed at the pres- 
ence of AF order. According to Eq.(|5^) V^= V^, = Yd 
and Eq.(^2|) can easily be diagonalised. One finally ob- 
tains the gap equations for the new eigenstates {k—+,-): 

di = -Y.[ilTbi:,)Vnm4Ff:, 

k' 

Ek^ie^ + dff^ (93) 
4 = (1 T 6£)^4 = (1 T 2a,/3,)^dl 
4 = -^(4±^4) 

These equations for pair states are formally 

identical to those for -00 , ^„ in Eq. (^ however with a 
different effective pair potential 

V±{k-k') = (1 T b^)^VD{k- k'){l T 6fe,)' (94) 

At the Bragg planes b, ~ 1, this implies that the pair 
potential V+ for d'^ is reduced to zero, similar as in the 

case K=0,u but for d^i the pair potential V_ is enhanced 
in this region of k- space, therefore AF order supports 
the sc pairing in the d^ state. 

As for the PM case in Sect. (V.A) we will now con- 
sider the simple situation of the FS sheet with cylindrical 
symmetry and small dispersion along c shown in the inset 
of Fig. H Then again the gap equations can be decom- 
posed into Fourier components which satisfy Eq. (^) as in 
the PM case. However the form factor am{kz) and the 
quasiparticle energy E(^,k2) are now different because 
they include the effect of AF order. One obtains from 
Eqs. (1^,11) 

K = 0: amikz) ^ g{kz)cos'^{mkz) 

E{^,k,) = [e + g{k,){Y,d:cosnk,f]i 

n 

K^u: arn{kz) = gikz) cos^{mkz) 

E{^,kz) ^ [e + g{kz){Y,dZsmnkzf]i (95) 
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K = ±: a,n{kz) = h^{k^)cos'^{mk:^) 

n 

Here we defined g{kz)—ai—{al^ — PlJ"^ and h±(k2 
(l=F2afe^/3fe^)- Using Eq.(pN) one can derive 



cos^ k. 



X2 



cos^ kz 
lT\l-9{k.)]^ 



(96) 



The functions g(kz), ^^(k^) with the parameter 
Aaf <C 1 given below accounts for the influence of AF 
order on the sc gap via the FS reconstruction caused by 



Ti^/ in Eq.(82). Specifically a finite Xaf will spit the 



degenerate doublet pair state (d„,d2,) of the PM case. 



D. Numerical solution of gap equations 

We now investigate some of the possible solutions of 
the gap equations quantitatively. It must be kept in mind 
that this work is based on a model theory which cannot 
make absolute predictions for the experimentally realised 
gap function. Its main purpose is a better understanding 
of the physics of a novel type of pairing mechanism lead- 
ing to the magnetic exciton mediated superconductivity. 

We first give estimates of the physical parame- 
ters involved in the present superconducting model 
for UPd2Al3 in addition to those already discussed in 
Sect. (Ill) for the magnetic exciton model itself. Accord- 
ing to Ref. [ Q the Fermi level is in the middle of a 
peak which we approximate by a square DOS of width 
W=2Ei^-=0.88 eV. Assuming one conduction electron per 
U site we have N(0)=e^^ = 0.0023 meV"^ for the model 
DOS. The smaU bandwidth W|| <W of the k^ disper- 
sion is given by W||=4t|| ~ 2Ej;^(A™''^/A^"'-l) where 
the last approximate relation is obtained assuming a free 
electron like dispersion _L c and relating the change of the 
FS cylinder cross section Ac in the PMBZ along c to the 
hopping t|| . The deviation of the ratio A""^^/A™*" from 1 
characterizes the amount of 'corrugation' of the FS cylin- 
der along c as sho"WH in Fig. ||, its value is obtained from 
dHvA experimentsta as A™''^/A™"' ~ 1.24 in the AFBZ. 
This leads to an estimate W||= 0.21eV or (W||/W)= 0.24 
consistent with the assumption made above. The inter- 
action He/ between localized 5f and conduction electrons 
with strength I=Q;Io(g-l) has two competing effects which 
are characterized by the dimensionless coupling constants 



A. 



(97) 



Whereas \sc characterizes the strength of the pair po- 
tential due to magnetic exciton exchange, \af is associ- 
ated with the effect of AF on the sc pair states. Since 



both are determined by the interaction constant I they 
are not independent and one has the relation 



1 (J)o W 

2 a 



[^(0)A]2A| 



(98) 



The quantities W||/W= 0.24 and [N(0)A]2= 0.12 are 
estimated within our model so that only one independent 
parameter Asc remains which will be fixed to achieve the 
proper Tc= 1.8 K. The solution of the gap equations 
Eq.([74|) for the AF case of Eq.(p5[) proceeds by iteration. 
As discussed in Sect. (V. A) the most favorable pair state, 
judging from the irreducible potentials in Fig. ^ is the 
odd parity doublet given by Eq. ([zol) . 




FIG. 8. Temperature dependence of gap function ampli- 
tudes d_(T) (full line) and d+(T) (dashed line) with AF or- 
der effect included [Xaf according to Eq.(^)). AAsc =4.3 
meV is chosen such that Tc=1.8 K. dashed-dotted line: AF 
moment effect suppressed {Xaf =0), with the same Asc, the 
transition temperature is then T°=1.27 K. inset: Schematic 
view of PM Fermi surface (FS) described by Eq.(|2|) with PM 
and AF Brillouin zones indicated. F and A' are the AFBZ 
and PMBZ zone centers respectively. States on the FS con- 
nected by the AF ordering wave vector Q are reconstructed 
by the AF order. 

The temperature variation of its sc amplitudes d+ (T) , 
d_(T) (i.f. the index m=l is suppressed) is shown in Fig. 
^ for two cases: (a) neglect of AF order, setting \ap =0 
arbitrarily then both representations arc accidentally de- 
generate (b) including the effect of AF with the proper 
\af according to Eq.(p8|): the degeneracy is hfted due to 
the finite Xaf and the nondegenerate d^ state describes 
the true sc phase in the AF background. It is prefered 
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as compared to due to its gap enhancement factor 
(1+bg) in Eg. (p3|) . The iteration procedure used for the 
calculation of ^^(k^) in Eq.(|7^) allows for the appearance 
of different Fourier components below Tck which corre- 
spond to V" irreducible potentials which are smaller than 
the maximum one with n=m. Their appearance would 
lead to a change of the profile of dni^z) as function of 
temperature. For all cases studied with the potentials of 
Fig. such a case never appears, therefore in the quasi- 
particle energies in Eq.(p5|) the summation over n can 
simply be replaced by the dominant term for n=m. 



(99) 



will appear. The k^- dependence of d^ and d^^^ 
{Xaf = 0) is illustrated in Fig. ^ (inset). Here AAsc= 
4.3 meV (Asc=0.72), which after Eq.(|98|) corresponds to 
A^F=0.018, is fixed to reproduce the physical value for 
Tp=1.8 K for d^. 



VI. A HYBRID MODEL INCLUDING SPIN 
FLUCTUATIONS 




I N(0) [meV] 



FIG. 9. Variation of transition temperature with coupling 
strength AAsc=I^N(0) for the (v,w) odd parity doublet states. 
With AF order effect suppressed {{Xaf =0) Tc's are degen- 
erate (dash- dotted line). With effect of Xaf (Eq.(^)) in- 
cluded a splitting appears (full and dashed lines). Therefore 
in the presence of AF only the nondegenerate d-{kz) state of 
Eq.(M) will be stable at any temperature below T" =1.8 K 
for AAi,c=4.3 meV. inset: wave number (k^-) dependence of 
the gap functions d±(k.z) and dv,w(kz) {Xaf =0) at T= 0.05 



K. The anomalies at the Bragg plane ^Qz 
effect of AF order. 



are due to the 



In Fig. ^ it is shown how the suppression of Tc and 
the splitting of d'^, d!l representations evolves as func- 
tion of Asc due to the effect of Xaf given in Eq. (|9^) . AF 
order renders the d^ state more favorable. Once d*i >0 

at Tc the split- off doublet component d^ with its lower 
transition temperature will be completely suppressed by 
the already finite condensation energy of the d^ state. 
Therefore in the AF background only a single nondegen- 
erate state 



It has been the aim of the previous analysis to inves- 
tigate a new type of pairing mechanism which is a con- 
sequence of the dual nature of 5f- electrons in UPd2Al3 . 
The model employed in Eq.(0) as an extreme case with- 
out any repulsion (lc=0) between 5f- conduction elec- 
trons of the type given in Eq.(|4|). Therefore only the 
new pairing mechanism via the exchange of magnetic ex- 
citons is possible which was then studied for its own sake. 
In a real 5f- compound one has to take into account also 
the qu asip article Coulomb repulsion (Ic >0) described 
by Eq. (|34| ) which leads to a conventional spin fluctuation 
type part for the total pairing potential. In this last sec- 
tion we therefore study a hybrid model where both contri- 
butions are simply added, neglecting any cross influence. 
We investigate the evolution of the pair potential and fa- 
vorable sc states as a function of the ratio p=(Ic/I) which 
controls the relative strength of these contributions. The 
pure magnetic exciton model studied before corresponds 
to p=0. The two contributions of Eqs.( ^4|j4^ ) to the total 
pair potential are determined by the RPA static conduc- 
tion electron susceptibility xio) and the exciton mode 
energy uj-{q) respectively. Due to Eqs. (|55| , |5^ ) they are 
related via 



(100) 



with ujQ ^ A' + Px, X = -wT.qX{c[)- Since the su- 
perexchange ^ ff{q) is unknown we restrict to the sim- 
plest possible case assuming ^ ff{q) —0, then Xoi't) is 
directly determined by the magnetic exciton dispersion 
^-{q) up to a constant wq which has to satisfy the con- 
straint 



I^Xo{q) = ujQ ~ uj-{q) > 



(101) 



Then, using Eqs.( |3qj5q ) the total pair potential of the 
hybrid model is given by 



2 ^R\q) 



2 u-{q) 



(102) 



1 
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where the RPA conduction electron susceptibihty of 
Eq.(^9|) is obtained as 



p^x{q) = 



1 - f [wo - 



(103) 
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FIG. 10. Effective pair potentials (absolute values) for the 
hybrid model for even (0) and odd (u,v,w) parity for p= Ic/I 
=0.124. inset: most important Fourier components as func- 
tion of model parameter p. With increasing p (increasing spin 
fluctuation contribution) the even parity singlet (k=0, m=l) 
becomes stable at p = 0.097 (first arrow) and then more fa- 
vorable than the odd parity u state (second arrow) but stays 
less favorable than the (v,w) odd parity doublet up to the 
maximum p= 0.124. 



For the hybrid model of Eq.(102) one can now calcu- 
late the (ix, qy- averaged potentials Y^ic^z) as in Eq.(|7l|) 
where the bar will now be suppressed. They are shown 
in Fig. ^ as f unct ion of qz for parameters P chosen 
such that Eq.(lOl) is fulfilled and stays paramag- 
netic, i.e. nonsingular. One observes that the singlet 
potential (V^) has now also a pronounced maximum at 
Q comparable to the potentials Yv,w for the odd parity 
doublet state. Note that for the hybrid model Vu ^ - 
Yv,w and therefore three potential functions are present 
in Fig. |l^ contrary to the situtation in the pure magnetic 
exciton model of Fig. The stability of the singlet sc 
state should now be increased as p increases. This can 
clearly be seen by plotting the relevant Fourier compo- 
nents V™ as function of increasing p, i.e. increasing the 
spin fluctuation contribution to the pairing. This is done 
for the unstable singlet state in Fig. ^ (dashed curve, 
m=l) versus the two most favorable odd parity states 



(up per arrows in Fig. ^ and the result is shown in Fig. 
|l^. For p > 0.09 one has -V^ >0 and the singlet state 
becomes possible. For even larger p > 0.12 it becomes 
more stable than the du{k) (m=2) odd parity state and 
then -Vs rapidly approaches the potentials Yy^w for the 
doublet state for th e lar gest value of p. At this value the 
model xio) in Eq.(lOl) becomes singular indicating an 
AF instability of the conduction electron system. Fo r the 
hybrid model with T ff{q)^0 and xio) given by Eq.(103) 
the odd parity state is therefore always the stable one, 
however the singlet state do(k2) becomes increasingly fa- 
vorable as the spin fluctuation part increases. In principle 
it is possible that in the gene ral case with the (unknown) 
J'//((j) included in Eq.(lOO) the hybrid model will pre- 
fer the singlet sc state do(fc). Of course, if we neglect 
the magnetic exciton part altogether and consider only 
the isotropic pure spin fluctuation model then accord- 
ing to Eq.(38) Ys ^Ya^-3Yu,v,w, and provided xio) has 
the maximum at q — Q, the singlet even parity sc state 
do(k.z)— djcosfc^ is the stable one. 



VII. SUMMARY AND OUTLOOK 

In this work a dual model of localised and itinerant 5f- 
electrons was studied as a generic model for superconduc- 
tivity in uranium- compounds and applied to UPd2Al3 . 
The dual nature of 5f electrons in this compound was 
concluded from experimental evidence, notably suscepti- 
bility and normal state Knight shift measurements. The 
localized electrons were assumed to be in a 5f^ configu- 
ration split by the CEF into singlet ground and excited 
states separated by an energy A=6 meV. The conduction 
electrons form band states whose main FS sheet can be 
described by a corrugated cylinder aligned along the c- 
axis. Due to RKKY and superexchange interactions the 
CEF excitation at A is broadened into a dispersive band 
of magnetic excitons which become soft at the Neel tem- 
perature Tn where an induced AF moment appears in 
the singlet ground state. These magnetic excitoB^pSSjhere 
found in inelastic neutron scattering experimentgaliJ and 
a theoretical explanation of their dispersion has been 
given within the AF singlet ground state model. |-J^ur- 
thermore a comparison with tunneling experimentst£l has 
suggested that the exchange of magnetic excit ns be- 
tween conduction electrons leads to a new type of pairing 
mechanism which is responsible for superconductivity in 
UPd2 AI3 . The consequences of this hypothesis have been 
analyzed in detail within the dual model. The new mag- 
netic exciton mediated pair potential was derived within 
a diagrammatic RPA type approach and compared to 
the results of the well known spin fiuctuation theory pair 
potential caused by anti- paramagnon exchange. It was 
found that the former manifestly violates rotational sym- 
metry in spin space contrary to the latter which is com- 
posed of spin singlet and triplet contributions due to ro- 
tational invariance. The new effective pair potential due 
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to magnetic exciton exchange has a lower symmetry with 
contributions from the even singlet, another odd parity 
nondegenerate and an odd parity doublet contribution. 
The projector representation defining the irreducible pair 
potentials has been derived and for the cylindrical FS 
sheet the effectively ID gap equations have been solved. 
The favoured representation is determined by the largest 
Fourier coefficient of the in-plane averaged pair poten- 
tial. In case of AF magnetism is excluded it was found 
that the new mechanism favors the odd parity doublet 
state which has equatorial node lines in the gap func- 
tion. The singlet state is generally not favored in the 
present weak coupling version of the magnetic exciton 
mechanism. The influence of singlet ground state AF on 
sc leads to important consequences: The doublet state is 
strongly split and the stable sc state is a nondegenerate 
odd parity state given by Eq. (|99|) . 

A hybrid model including spin fluctuations has also 
been considered. In the simplified case that superex- 
change between localized 5f electrons can be neglected 
the spin fluctation pair potential can also be obtained 
from the magnetic exciton dispersion. As expected its 
admixture leads to a stabilization of the singlet pair state 
with node lines at ±^Q. However within the simplified 
model the above mentioned odd parity state remains the 
sc ground state. In a more general hybrid model which 
does not have a direct connection between the two con- 
tributions to the pair potential an even singlet sc ground 
state is possible which is favored for the pure spin fluc- 
tuation model of the pair potential. In this work the 
thermodynamic signatures of the proposed and possible 
gap functions have not been discussed. Their detailed 
investigation is an important aspect foiJulure work. 

Thermal conductivity measurementsEZrtj point to the 
existence of node lines in the gap in a plane parallel to 
the hexagonal ab plane but sofar do not allow to give 
their positipH ck^ along the c axis. Recently it has been 
proposedcal23 that the field-angle dependence of thermal 
conductivity gives a unique possibility to detpcminc the 
position of the node lines. It has been shownL3 that the 
angle dependence should be distinctly different for an 
odd parity gap function sinfc^ with k^=0 node line as in 
Eq. ( |79[ ) and the even parity gap function cos with node 
line at kz=-| as in Eq.(pO|). It is generally hdieiffid that 
the observation of a Knight shift below tJ3'EJo points 
to an singlet even parity gap function as in Eq. (^0|) which 
is not favored by the present form of the magnetic exci- 
ton model. However the Knight shift has also large con- 
tribution from the AF ordered localised moments and 
the quantitative extraction of the itinerant part is not 
unique. Therefore one has to await the field- angle de- 
pendent thermal conductivity measurements to solve the 
problem of the true node line position and parity of the 
order parameter. 

The theory developed in this paper illustrates the ba- 
sic mechanism of Cooper pair formation via the magnetic 
exciton exchange in U-compounds with dual (localised 



and itinerant) 5f electrons. In its context differences to 
the conventional spin fluctuation mechanism and general 
properties of the nonretarded pair potentials were inves- 
tigated. However one cannot expect fully quantitative 
predictions for UPd2Al3 from such a simplified weak cou- 
pling model. Ultimately a strong coupling theory with 
a retarded magnetic exciton mediated pair potential is 
necessary for a realistic application of this new pairing 
mechanism to UPd2Al3 and a theoretical approach to the 
question of stability of even or odd parity gap function. 
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APPENDIX A: 



APPENDIX B: 



In this appendix we give the basic geometry of the con- 
ventional unit cells in the direct (Wigner Seitz cell: WS) 
and reciprocal (Brillouin zone: BZ) hexagonal lattice for 
the discussion of the magnetic exciton dispersion. The 
BZ is rotated by 30° with respect to the WS cell. The 
relevant primitive lattice vectors of the direct and recip- 
rocal lattices are given, respectively, by 



a( 
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2' 



x/3 



,0) 



b* 
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a(2,— ,0) 
c(0,0,l) 



— 0,0,1 
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(Al) 
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Any wave vector can be written in this basis as q—ha* + 
kb*+ Ic* . In reduced units of ^ and ^, 9 is then given by 



PMBZ 
AFBZ : 



q = 2Tr{h + k,^{h~k),l) 



<f = 27r{h + k,^{h^k),^l') 



(A3) 

For the qz components the relation I' = 21 — 1 or I — 
^{I' + l) holds, i.e. the AFBZ is obtained from the PMBZ 

through the shift by an AF wave vector Q = (0, 0, tt) and 
scaling by a factor i. The corresponding wave vectors 

q ' = q — Q of the symmetry points of the hexagonal 
AFBZ are then given by 

r : 7r(0,0,0) 
M:7r(l,-i=,0) 

if :7r(|o,0) 



A:7r(0,0,-) 
A' : 7r(0,0, 1) 



L:7r(l, 



1 1 

V3'2 
4 1 



In the AFBZ the symmetry point A' is equivalent to F 
since it is shifted by a reciprocal lattice vector Q. 



The single- ion susceptibility tensor which is a 

necessary ingredient for the calculation of the magnetic 
exciton dispersion from Eq. (^8|) is obtained from the gen- 
eral linear response expression for two operators A, B in 
a system of two levels e± corresponding to CEF singlet 
eigenstates |±) in the molecular field: 



uba{lo) ^ P{T)[ 
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P{T) = 2(5) = tanh|A' 
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BA 
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6+ — e_ is the CEF excitation 
: ± refer to the mf- eigenstates 



Here {S) = (S, 
energy of Eq.(|l2|) 
in the ordered phase (Eq.(|llj)) and A=A-{A), B=A-{B). 
For A,B= Jq (a =x,y,z) the form of the 3x3 dipolar 
cartesian susceptibility tensor ??(w) is too cumbersome 
for solving the secular equation for the exciton modes. 
This may be remedied by performing a real space rotation 
(x,y,z)^(x',y,z') around the y- axis by an angle 9r such 
that it exactly compensates the effect of the rotation in 
state space which transforms to the eigenstates |±) of the 
mf- Hamiltonian Eq. (|ll|) . In this basis one has 
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with a=cos2^^=[l-|-7'^] i ^h—sm2d — j'a. Performing 
the counter- rotation leads to (u^ = cos^r, v,. = sin^^) 

Jx ^rJx ^rJz: Jy — Jy i — ^rJx ^rJz (-E^^) 

defining (j> = 29 this can be written 
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sm(t/ 
— cos( 



By choosing 0^ = —4> = —20 the two transformations 
compensate and Tx—^x, ■i\=3y, J'z=Jz- Therefore in 
the rotated (x',y',z')- coordinate system the susceptibility 
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tensor of the AF ordered two level system is identical to 
the same tensor for the paramagnetic case in the original 
(x,y,z) hexagonal coordinate system. It has only nonzero 
elements for x,y and is given by 



u'icu) 



A' 



A' 
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A' -iuj 

ILU A' 



This form will be used in the secular equation (|20|) 



(B5) 



^ M. Sigrist and K. Ueda, Rev. Mod. Phys. 63, 239 (1991) 

2 R. Fisher, S. Kim, B. Woodfield, N. Phillips, L. TaiUefer, K. 
Hasselbach, J. Flouquet, A. Giorgi and J.L. Smith, Phys. 
Rev. Lett. 62, 1411 (1989) 

3 G. Bruls, D. Weber, B. Wolf, P. Thalmeier, B. Liithi, A. 
de Visser and A. A. Menovsky, Phys. Rev. Lett. 65, 2294 
(1990) 

"^F. Kromer, R. Helfrich, M. Lang, F. Steglich, C. Lang- 
hammer, A. Bach, T. Michels, J.S. Kim and G.R. Stewart, 
Phys. Rev. Lett. 81, 4476 (1998) 

^ H. Amitsuka, M. Sato, N. Metoki, M. Yokoyama, K. Kuwa- 
hara, T. Sakakibara, H. Morimoto, S. Kawarazaki, Y. 
Miyako and J. A. Mydosh Phys. Rev. Lett. 83, 5114 (1999) 
Ch. Geibel, C. Schank, S. Thies, H. Kitazawa, C. D. Bredl, 
A. Bohm, M. Ran, A. Grauel, R. Caspary, R. Helfrich, U. 
Ahlheim, G.Weber and F. Steglich, Z. Phys. B 84, 1 (1991) 
R. Grauel, A. Bohm, H. Fischer, C. Geibel, R. Kohler, R. 
Modler, C. Schank, F. Steglich, G. Weber, T. Komatsubara 
and N. Sato, Phys. Rev. B 46, 5818 (1992) 

* R. Feyerherm, A. Amato, F.N. Gygax, A. Schenck, C. 
Geibel, F. Steglich, N. Sato and T. Komatsubara, Phys. 
Rev. Lett. 73, 1849 (1994) 

® T.E. Mason and G. Aeppli, Matematisk-fysiske Med- 

delelser 45, 231 (1997) 
^° N. Bernhoeft, N. Sato, B. RoessU, N. Aso, A. Hiess, G.H. 
Lander, Y. Endoh and T. Komatsubara, Phys. Rev. Lett. 
81, 4244 (1998) 

" N N. Sato, N. Aso, G.H. Lander, B. RoessU, T. Ko- 

matubara and Y. Endoh, J. Phys. Soc. Jpn. 66, 1884 (1997) 

and 66, 2981 (1997) 
^2 N. Bernhoeft, Eur. Phys. J. B 13, 685 (2000) 
" M. Jourdan M. Huth, H. Adrian, Nature 398, 47 (1999) 
" N.K. Sato, N. Aso, K. Miyake, R. Shiina, P. Thalmeier, 

G. Varelogiannis, C. Geibel, F. Steglich, P. Fulde and T. 

Komatsubara, Nature 410 340 (2001) 

R. Shiina, private communication 

M. Jourdan M. Huth, H. Adrian, Nature 398, 47 (1999) 
^■^N.K. Sato, N. Aso, K. Miyake, R. Shiina, P. Thalmeier, 
G. Varelogiannis, C. Geibel, F. Steglich, P. Fulde and T. 
Komatsubara, Nature 410 340 (2001) 

K. Knopfle, A. Mavromaras, L.M. Sandraatski, J. Kiibler, 
J. Phys. Cond. Matter 8, 901 (1996) 



Y. Inada, H.Yamagami, Y. Haga, K. Sakurai, Y. Tokiwa, 
T. Honma, E. Yamamoto, Y. Onuki and T. Yanagisawa J. 
Phys. Soc. Jpn. 68, 3643 (1999) 

J. Jensen and A.R. Mackintosh, 'Rare Earth Magnetism, 
Structures and Excitations' Clarendon Press, Oxford 1991 
A. Krimmel, P. Fischer, B. Roessh, H. Maletta, C. Geibel, 
C. Schank, A. Grauel, A. Loidl and F. Steglich, Z. Phys. B 
86, 161 (1986) 

S.Nakajima, Progress of Theoretical Physics 50, 1101 
(1973) 

K. Miyake, S. Schmitt-Rink and CM. Varma Phys. Rev. 
B 34, 6554 (1986) 

A. L Morosov, JETP 83 ,1048 (1996) 

A. Amici P. Thalmeier and P. Fulde, Phys. Rev. Lett. 84, 
1800 (2000) 

A. Amici, P. Thalmeier and P. Fulde, Physica C 317-318, 
471 (1999) 

M.Chiao, B. Fussier, B.EUman, L.Taillefer, Physica B 230- 
232, 370 (1997) 

M. Hiroi, M. Sera, N. Kobayashi, Y. Haga, E. Yamamoto 
and Y. Onuki, J. Phys. Soc. Jpn. 66, 1595 (1997) 
Y. Sun and K. Maki, Euro phys. Lett. 32. 355 ( 1995) 
P. Thalmeier and K. 



Maki, 



cond-mat/0109351 



H. Tou, Y. Kitaoka, K. Asayama, C. Geibel, C. Schank and 
F. Steglich, J. Phys. Soc. Jpn. 64, 725 (1995) 



32 



K. Matsuda, Y. 
15223 (1997) 



Kohori and T. Kohara, Phys. Rev. B 55, 



21 



